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O ABSTRACT 

<N 

*^ Context It has been shown that by means of different physical mechanisms the expansion of H n regions can trigger the formation of new stars of 
jp* all masses. This process may be important to the formation of massive stars but has never been quantified in the Galaxy. 
^H Aims. We use Herschel-PACS and -SPIRE images from the Herschel Infrared survey of the Galactic plane, Hi-GAL, to perform this study. 
^J Methods. We combine the Spitzer-GUMPSE and -MIPSGAL, radio-continuum and sub-millimeter surveys such as ATLASGAL with Hi-GAL 
^«. to study Young Stellar Objects (YSOs) observed towards Galactic Hn regions. We select a representative Hn region, N49, located in the field 
centered on /=30° observed as part of the Hi-GAL Science Demonstration Phase, to demonstrate the importance Hi-GAL will have to this field of 
research. 
^_^ Results. Hi-GAL PACS and SPIRE images reveal a new population of embedded young stars, coincident with bright ATLASGAL condensations. 
_^f4 The Hi-GAL images also allow us, for the first time, to constrain the physical properties of the newly formed stars by means of fits to their spectral 
^_^ energy distribution. Massive young stellar objects are observed at the borders of the N49 region and represent second generation massive stars 
\^ whose formation has been triggered by the expansion of the ionized region. 

Conclusions. The first Hi-GAL images obtained using PACS and SPIRE have demonstrated the capability to investigate star formation triggered 
by H ii regions. With radio, submillimeter, and shorter wavelength infrared data from other surveys, the Hi-GAL images reveal young massive 
star- forming clumps surrounding the perimeter of the N49 Hn generated bubble. Hi-GAL enables us to detect a population of young stars at 
q different evolutionary stages, cold condensations only being detected in the SPIRE wavelength range. The far IR coverage of Hi-GAL strongly 
S^ constrains the physical properties of the YSOs. The large and unbiased spatial coverage of this survey offers us a unique opportunity to lead, for 
£^J the first time, a global study of star formation triggered by H n regions in our Galaxy. 
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^ 1. Introduction their borders. This indicates that triggering is important in the 

._ creation of massive stars and that hot photodissociation regions 

Ionized (H n) regions are known to trigger the formation of stars (pDRs) are a good place t0 study the earliest phases of massive . 

^i by means of various physical mechanisms (see Elmegreen ll998l star f orma ti n 
. and Deharveng et al. |2005| for a review). Several Hn regions 

IT} have been studied individually in the context of triggered star Our long-term aim is to use Hi-GAL (Molinari et al. l20T0T >. 

O formation, focusing on the associated neutral material and the combined with existing infrared, submillimeter, and radio sur- 
young stellar population (Zavagno et al. 2006 ; Deharveng et V eys, to study the influence of H n regions on triggering the fer- 
al. 2009; Pomares et al. 2009; Bieging et al. |2009) . These stud- mation of stars in our Galaxy. Hi-GAL's extended wavelength 
j>! ies have shown that the expansion of H n regions can trigger the coverage towards the far-infrared and its unprecedented sensitiv- 

• »-j formation of new stars of all masses. The Sp itzer-G UMPSE ity offer a unique opportunity to detect an embedded population 

/\ survey of the Galactic plane (Benjam in et al. I2003I I detected f young sources that are not detected at shorter wavelengths. 

H nearly 600 bubbles (Churchwell et al. I2006l >. Deharveng et al. This allows us to observe intermediate and high-mass YSOs over 

020 lOl l selected a series of 102 ionized bubbles and studied the the complete range of evolutionary stages. The unprecedented 

star formation in their su rround ings using Spteer-GLIMPSE and resolution of Hi-GAL also offers the opportunity to accurately 

MIPS GAL (Carey et al. |20 Q9J, ra dio (MAGPIS; Helfand et al. characterize the physical nature of the sources by means of a 

120061 and V GPS; Stil et al. !20061 >, and ATLASGAL (Schuller detailed fit to their spectral energy disttibution (SED). 
et al.[2009) data. They show that 86% of these bubbles enclose 

H ii regions, and that more than 20% of 64 bubbles (for which In this Letter we study the bubble-shaped ionized region, 
the ATLASGAL angular resolution is sufficient to resolve the N49, from the Churchwell et al. (2006) catalogue to illustrate the 
spatial distribution of cold dust) show massive star formation on purpos e of o ur project. The N49 bubble was s tudied by Watson 
et al. (2008) and by Deharveng et al. (120101 ). However, these 



* Herschel is an ESA space observatory with science instruments pro- studies had no information in the 70-500 fim range. This infor- 
vided by European-led Principal Investigator consortia and with impor- mation, obtained with the PACS and SPIRE data presented here, 
tant participation from NASA. allows us to discuss the star formation of this region in detail. 



A. Zavagno et al.: Star formation triggered by H n regions in our Galaxy 




Fig. 1. N49: a) Three-color composite image with 8 /mi Spitze r-GUMPSE (blue), 24 /mi MIPSGAL (green), 20 cm MAGPIS (red). 
ATLASGAL 870 /mi contours are superimposed on all images (levels are 0.01, 0.1, 0.22, 0.31, 0.5, and 2 Jy/beam). ATLASGAL 
condensations 1 to 4 are identified. We have identified the two massive stage I YSOs (YSO#l and YSO#3) discussed by Watson et 
al. J2008). The ultracompact Hn region is also identified; b) PACS 70 /mi image; c) PACS 160 /mi image; d) SPIRE 250/rni image; 
e) SPIRE 350 /mi image; f) Two-color composite image with PACS 70//m (blue) and SPIRE 350/mi image (red). The field is 5' x 
5', centered on /=28.83°, £=-0.23° 



2. Hi-GAL results on N49 

We use PACS and SPIRE images obtained in parallel mode at 70, 
160, 250, 350 and 500 /mi (at resolutions of 5", 11.4", 17.9", 
25"and 35.7", respectively) to determine the influence of Hn 
regions on their surroundings using the N49 Hn bubble as an 
example. Our general goal is to quantify the efficiency of H n re- 
gions on triggering star formation around the perimeter of their 
radiation-driven expanding bubble. Figure [JJ shows images of 
this region at different wavelengths, including the new PACS and 
SPIRE images. ATLASGAL 870 //m contours are superimposed 
on all images. This H n region is centered on Z=28.83°, b=-0.23° 
and located at 5.5 kpc (Anderson et al. 2009). It has a diameter 
of 4 pc. A compact Hn region lies on its border (see Fig. [Tk). 
ATLASGAL 870 /mi emission shows a fragmented shell of col- 
lected material surrounding the ionized region (green contours 
on Fig. [Til. Three condensations are seen surrounding the ion- 
ized region and a secondary peak is detected farther away. These 
are referred to as condensations 1 through 4, respectively, and 
are indicated in Fig. [Tk. Using ATLASGAL data, Deharveng et 
al. (120 1 Oi l derived a mass for the entire shell of 4250 M Q , and 
the mass of the brightest condensation in this shell (condensa- 
tion 1 in Fig. [TV) of 2350 Mq. The NH3 velocity measurements 
show that the condensations 1, 2, and 3 are associated with the 
ionized region (Wienen et al., in preparation). No measurement 



exists for condensation 4 but 13 CO (1-0) data from the Galactic 
Ring survey (Jackson et al. 2006) show that it is also associ- 
ated with N49. Two young stellar objects (identified YSO#l and 
YSO#3 in Fig. [Tk) are seen as bright 24 /mi sources towards 
condensation 1, the brightest at 870 /mi. These two YSOs are 
massive and in an early evolutionary stage (Watson et al. 2008). 
Using high resolution 6.7 GHz methanol maser measurements, 
Cyganowski et al. (2009) suggest that a rotating disk is associ- 
ated with YSO#3. An ultracompact H 11 region nearby YSO#3 is 
shown in Fig. [Tk as a bright, but more diffuse, source at 24 /mi. 
According to Deharveng et al. (120101 ). this region may have been 
ionized by a B0V star. 

Figure [TV) shows the PACS 70 /mi image with superimposed 
ATLASGAL 870/mi contours. We see that YSO#3 dominates 
the emission at longer wavelengths. The ultracompact Hn re- 
gion, bright at 24 /mi on the MIPSGAL image, is not bright in 
the 70-500 //m range. A 70 /mi source is detected towards con- 
densation 2 and 3. These sources are not detected at 24 /mi. 

Figure[JJl shows the SPIRE 250 /mi image with ATLASGAL 
870 /mi contours superimposed. The sources previously detected 
at 70 /mi towards condensations 2 and 3 are seen. Weak emission 
is observed towards condensation 4. This emission was not de- 
tected at 70 //m and 160 /mi. 

The PACS images clearly show the hot PDR that surrounds 
the H 11 region. The 70 //m and the 8 /mi emission of the PDR are 
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very similar. The only difference comes from the PACS 70 /mi 
source seen towards condensation 2 and observed as a dark spot 
on the 8 /mi image. This indicates that condensation 2 and the 
PACS IR source are probably located in front of the PDR. The 
situation is different for condensation 3, for which no clear ab- 
sorption is seen at 8 /mi suggesting that the condensation is lo- 
cated behind the PDR. The two color-composite image shown 
in Fig. []}, PACS 70 /«n (blue) and SPIRE 350 /an (red), indi- 
cate that parts of the region where star formation is occurring 
are denser or colder or both, as indicated by the red emission. 

The brightest IR source observed around N49, probably cor- 
responding to YSO#3, is located along an IR dark cloud seen 
on the 8 /mi GLIMPSE image. Since blobs are observed along 
this filament on larger scales, the origin of the formation of this 
bright IR source is questionable. Was it produced by the frag- 
mentation and collapse of the shell surrounding N49 by means 
of the collect and collapse process (Elmegreen & Lada |1977] l, or 
by the compression of a pre-existing blob by the N49 (and the ul- 
tracompact) Hn region(s)? Morphological considerations favor 
the first option as a pre-existing clump would have distorted the 
ionization front, which is not what we observe here. We discuss 
the star formation observed around N49 in the framework of the 
collect and collapse model in Sect. 3. 

2. 1 . The spectral energy distribution fitting 



We used the on-line SED fitting tool of Robitaille et al. (12007b 
to derive the physical properties of the sources observed towards 
ATLASGAL condensations bordering N49. Fluxes were mea- 
sured at 70, 160, 250, 350, 500, and 870 /mi by rebinning all the 
images to the PACS 70 //m resolution, using a square aperture of 
100", 50", 70", and 40" for condensations 1, 2, 3, and 4, respec- 
tively (2.67, 1.33, 1.87, and 1 pc for the N49 distance of 5.5 kpc), 
and subtracting emission from a local background. 

Watson et al. (2008) identified seven YSOs towards N49. 
Their YSO#5 and #6, observed near condensation 4, are prob- 
ably not related because they are not observed at longer wave- 
lengths on the Herschel images. Watson et al. (2008) did not 
detect the red sources observed with PACS and SPIRE to- 
wards the ATLASGAL condensations 2 and 3. As seen in their 
Fig. 17, their SED fits are not constrained at wavelengths beyond 
24 /mi due to a lack of data. The Hi-GAL coverage between 70 
and 500 /mi, together with the ATLASGAL 870 /<m emission, 
strongly constrain the SED of the observed sources, allowing us 
to derive their physical properties. 

In Table 1, we present our results of the SED fitting for all 
sources, taking into account fits with^ 2 ~xi t P er datapoint < 3. 



Table 1. Results of the SED fitting 



(M Q ) 



M, 



envelope 



■'"envelope 

(M Q ) 



L 
(10 3 Lq) 



Cond. 1 20-30 
Cond. 2 8-10 
Cond. 3 8-15 



8-10 

4 
2-5 



5000 



2000-5000 



50 
1 

2-5 



Towards condensation 1, a high luminosity, high mass source 
is found that has a high accretion rate. The values given in 
Table 1 are for the entire condensation 1, including YSO#l and 
YSO#3. The 70/zm emission is dominated by YSO#3. This is an 
extended green object and is associated with a methanol maser 




1 000 



Fig. 2. Results of the fit to the fluxes measured towards conden- 
sation 1 in N49 using the on-line SED fitting tool of Robitaille 
et al. (2007). The best-fit (dark solid line) is strongly constrained 
by the Hi-GAL wide wavelength coverage. The grey lines show 
all the range of models. The 24 /mi flux is taken from Watson et 
al. (2008). 



(Cyganowski et al. 2008, 2009), which implies a massive YSO 
with an age younger than 3.5xl0 4 years (Breen et al. 120101 1. 
However, no emission is seen in its direction at 20 cm on the 
MAGPIS image. The high accretion rate derived for condensa- 
tion 1 (see Table 1) may prevent the development of the ionized 
region (Churchwell 2002). The emission may also be optically 
thick at 20 cm. The results of the fit for condensation 1 are shown 
in Fig. [2] 

The sources observed towards condensations 2 and 3 are 
also massive objects with high accretion rates. The emission 
that is barely detected at 250 //m towards condensation 4 can- 
not be fit using the Robitaille et al. SED fitting tool because no 
data at shorter wavelengths exist to constrain the fit. This ob- 
ject is probably a cold core with no internal source of heating. 
However, it is interesting to detect this emission farther away 
from the ionization front as it may represent an earlier evolution- 
ary stage in the star formation process. The 870/mi ATLASGAL 
peak flux of 0.26 Jy/beam indicates a column density of N(H2)= 
6.8xl0 21 cm" 2 (2.2xl0 22 cm" 2 ) for a temperature of 20 K 
(10 K). This represents a low density condensation compared 
to the values derived by Deharveng et al. (2010) around Galactic 
bubbles where star formation is clearly observed. 

In all cases, the age of the source is not well constrained by 
the fitting procedure. Other indicators, such as chemical clocks, 
outflows, and masers, are needed to more tighly constrain the 
ages. 

3. Discussion 



Following the analytical approach of Whitworth et al. (1994 see 
their Sect. 5), we discuss the observed star formation around 
N49 in the context of the collect and collapse process. The 
present radius of N49 (2 pc), its estimated age (0.5-1 Myr), and 
the spectral type of its ionizing star (05 V, Watson et al. 2008 
Everett & Churchwell 2010) allow us to derive a mean value for 
the initial density of the surrounding medium, n =* 10 4 cm -3 . 
Assuming this density and following the formula given by 
Whitworth et al. (1994), we can estimate the time at which the 
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fragmentation begins, ff rag ment -0.5 Myr and the radius and col- 
umn density of the shell at this time, Rf mgmen i -1.55 pc and 
-^fragment -1.6xl0 22 cnT 2 , which are similar to the values we 
derived towards condensation 4. All these values strongly de- 
pend on the sound speed in the shocked gas a, (which is between 
0.2 and 0.6 km/s) and n , the initial density of the medium. The 
values we give here are derived using a,=0.2 km/s and no=10 4 
cm -3 . If no is lower, we obtain a f fragment higher than the esti- 
mated age for N49, which is not possible because the fragments 
are already seen. The same happens if we increase a,. However, 
the first set of values that provide acceptable results for the start 
time of the fragmentation and radius (and column density) of 
the shell at that time does not provide acceptable results for the 
mass of the fragments, which are too low (around 7.5 M ) com- 
pared to the masses derived by Deharveng et al. ( 120101 ) using 
the ATLASGAL measurements (of about 2500 M Q , 240 M , and 
355 M for condensations 1, 2 and 3, respectively). As discussed 
below, the wind dynamics that exists towards N49 can mod- 
ify this simplified approach. A finer age determination for the 
sources observed towards the ATLASGAL condensations may 
help in more tightly constraining this model. 

The Hi-GAL images and results of the SED fits show that 
five massive YSOs (YSO#l and #3, the UC Hn region, and 
IR sources towards condensations 2 and 3) are formed (or are 
in the process of forming) on the border of N49. This region 
is ionized by an 05V star and stellar winds are believed to 
play a crucial role in the formation of this bubble. Everett & 
Churchwell (120101) studied in detail the case of N49 as a wind- 
blown bubble. They conclude that the survival of dust within 
the ionized gas in the presence of winds requires a high den- 
sity in the surrounding medium. This high density together with 
greater dynamical motion caused by stellar winds might have 
increased the efficiency of triggered star formation, in terms of 
number and mass of the second-generation stars. Other dusty 
wind-blown bubbles will be observed in Hi-GAL and the unbi- 
ased access to the YSOs population will help us to understand 
whether winds play a crucial role in increasing the efficiency of 
triggered star formation. 



4. Conclusions 

We have presented the first PACS and SPIRE images from Hi- 
GAL of the bubble-shaped Galactic H n region N49. This region 
is used as an illustration of the study dedicated to the star forma- 
tion triggered by Galactic H n regions that we plan to lead for the 
whole survey, combining the Hi-GAL results with other infrared 
and radio surveys of the Galactic plane. We have shown that: 

- The Hi-GAL SPIRE and PACS images allow us to study the 
distribution of young sources towards N49. The far-IR fluxes 
have been measured and strongly constrain the spectral en- 
ergy distribution of these sources. This allows us to charac- 
terize their properties. 

- The PACS images reveal the existence of red young stel- 
lar objects towards two ATLASGAL condensations, sources 
that had not been previously detected at shorter wavelengths. 
The ultracompact Hn region is not clearly seen in the 
Herschel range. The bright YSO#3 Watson et al. observed to 
be coincident with condensation 1 dominates the emission at 
longer wavelengths. 

- SED fits for the 3 sources detected by Herschel towards 
millimeter condensations using the Robitaille et al. (2007) 
model show that these sources are young and massive. 
However, their age has not been constrained and other in- 



dicators are needed to refine the discussion of star formation 
history in this region. 
- Five massive stars are forming in the N49 PDR. The high 
star formation efficiency in N49 may be due to the presence 
of winds from the first generation massive star. 

The study of N49 using a multiwavelength approach shows that 
Hi-GAL enables measurements for the crucial far-IR range to 
be made that are essential to infer the properties of the YSOs. 
Hi-GAL clearly provides important insight into star formation 
triggered by expanding Hn regions. Seventy-six Hn regions of 
all shapes are detected in the /=30° field and a higher density 
of Hi-GAL sources is clearly observed towards these regions, 
indicating that star formation triggered by H n regions may be 
an important process. 

Acknowledgements. Part of this work was supported by the ANR {Agence 
Nationale pour la Recherche) project "PROBeS", number ANR-08-BLAN-0241 
(LA). 



References 

Anderson, L., Bania, T.M. 2009, ApJ, 690, 706 

Benjamin, R.A., Churchwell, E., Babler, B.L., et al. 2003, PASP, 115, 953 

Bieging, J.H., Peters, W.L., Vilaro, B.V., Schlottman, K., Kulesa, C. 2009, AT, 

138, 975 
Breen, S.L., Ellingsen, S.P., Caswell, J.L., Lewis, B.E. 2010, MNRAS, 401, 2219 
Carey, S.J., Noriega-Crespo, A., Mizuno, D.R., Shenoy, S., Paladini, R., 

Kraemer, K.E., Price, S.D., Flagey, N., Ryan, E., Ingalls, J.G. et al. 2009, 

PASP, 121,76 
Cyganowski, C.J., Whitney, B.A., Holden, E., Braden, E., Brogan, C.L., 

Churchwell, E., Indebetouw, R., Watson, D.F., Babler, B.L., Benjamin, R. 

etal. 2008, AJ, 136, 2391 
Cyganowski, C.J., Brogan, C.L., Hunter, T.R., Churchwell, E. 2009, ApJ, 702, 

1615 
Churchwell, E. 2002 ARA&A, 40, 27 

Churchwell, E., Povich, M. S., Allen, D., et al. 2006, ApJ, 649, 759 
Deharveng, L., Zavagno, A., Caplan, J. 2005, A&A, 433, 565 
Deharveng, L. et al. 2009, A&A, 496, 177 

Deharveng, L., Schuller, E, Anderson, L. et al. 2010, A&A submitted 
Elmegreen, B.G., & Lada, C.J. 1977, ApJ 214, 725 
Everett, J.E., Churchwell, E. 2010, ApJ, in press 
Elmegreen, B.G. 1998, in ASP Conf. Ser. 148, ed. C.E. Woodward, J.M. Schull, 

& HA. Tronson, 150 
Helfand, D.J., Becker, R.H., White, R.L., Fallon, A., Tuttle, S. 2006, AJ, 131, 

2525 
Jackson, J.M. et al. 2006, ApJS, 163, 145 
Molinari, S. et al. 2010, A&A this volume 
Pomares, M., Zavagno, A. et al. 2009, A&A, 494, 987 

Robitaille, T., Whitney, B.A., Indebetouw, R., Wood, K. 2007, ApJS, 169, 328 
Schuller, E, Menten, K.M., Contreras, Y. et al. 2009, A&A, 504, 415 
Stil, J.M., Taylor, A.R., Dickey, J.M. et al. 2006, AJ, 132, 1 158 
Watson, C, Povich, M.S., Churchwell, E.B., et al. 2008, ApJ, 681, 1341 
Whitworth, A.P., Bhattal, A.S., Chapman, M.J., Turner, D., Turner, JA. 1994, 

MNRAS, 268, 291 
Zavagno, A., et al. 2006, A&A, 446, 171 



1 Laboratoire d'Astrophysique de Marseille (UMR 6110 CNRS & 
Universite de Provence), 38 rue F. Joliot-Curie, 13388 Marseille 
Cedex 13, France e-mail: annie . zavagno@oamp . fr 

2 Astrophysics Research Institute, Liverpool John Moores University, 
Twelve Quays House, Egerton Wharf, Birkenhead CH41 1LD 

3 Center for Astrophysics and Space Astronomy Department of 
Astrophysical and Planetary Sciences 389 UCB Boulder Colorado 

4 Spitzer Postdoctoral Fellow, Harvard-Smithsonian Center for 
Astrophysics, 60 Garden Street, Cambridge, MA, 02138, USA 

5 School of Physics, University of Exeter, Stacker Road, Exeter, EX4 
4QL, UK 

6 Max-Planck Institut fur Radioastronomie, Auf dem Hiigel 69, D- 
53121 Bonn, Germany 

7 ESO, Karl Schwarzschild-Strasse 2, D-85748 Garching bei 
Miinchen, Germany 



A. Zavagno et al.: Star formation triggered by H n regions in our Galaxy 



8 NASA Herschel Science Center, IPAC, Caltech, Pasadena, CA 
91125, USA 

9 INAF-IFSI, Fosso del Cavaliere 100, 00133 Roma, Italy 

10 School of Physics and Astronomy, Cardiff University, Queens 
Buildings, The Parade, Cardiff, CF24 3AA, UK 

11 Dipartimento di Fisica, Universit di Roma 2 "Tor Vergata", Rome, 
Italy 

12 ASI Science Data Center, 1-00044 Frascati (Rome), Italy 



